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Abstract. We study the global existence and decay to spherical equilibrium 
of Hele-Shaw flows with surface tension. We prove that without injection of 
fluid, perturbations of the sphere decay to zero exponentially fast. On the 
other hand, with a time-dependent rate of fluid injection into the Hele-Shaw 
cell, the distance from the moving boundary to an expanding sphere (with 
time-dependent radius) also decays to zero but with an algebraic rate, which 
depends on the injection rate of the fluid. 
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1. Introduction 

1.1. The problem statement. We establish decay estimates for solutions of the 
Hele-Shaw equations with surface tension on the free-boundary, and with injection 
of fluid into the cell. With the timc-dcpcndcnt fluid domain denoted by £l(t), an 
open subset of M 2 with boundary T(t), and for time t £ [0, T], the two-dimensional 
Hele-Shaw equations are given by 

Ap = -nd in Q(t) , (la) 

p = R on T(t), (lb) 

v(r(*)) = -j| on r(t), (ic) 

O(0) = f2, (Id) 

where S is the Dirac delta function at the origin, p(x, t) denotes the fluid pressure, 
/i = fx(t) denotes the rate of injection of fluid if fj, > (or suction if fj, < 0), and 
H is the (mean) curvature of the evolving free-boundary T(t). We use V(T(t)) to 
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denote the normal velocity of the moving free-boundary P(i), and we let n denote 
the outward-pointing unit normal on T(t). 

When the injection rate p ^ 0, the volume of the fluid domain |Sl(i)| can be 
computed as 



t 

2 



|fi(t)| =tt+ / p( S )ds = 7rp(t) 2 , (2) 



where pit) = \ 1+ [ ^^-ds > is the radius of a ball centered at the origin. 
V Jo ™ 

In fact, if the initial domain SI = B\ = £?(0,1), then the solution to (JTJ) is given 
explicitly by (l(t) = B(0, p(t)), with pressure function 

. . 1 pit) , led 

We will show that under certain growth conditions on the injection rate p(t), if 

(1) the initial domain SI is sufficiently close to the unit ball B\ = B(0, 1) in M 2 
with |0| = = tt; 

(2) the center of mass of SI is sufficiently close to the origin (which is the point 
of injection), 

then (l(i) converges to B p = B(0, p(t)) as t — > oo. The precise statement of our 
result is given below in Theorem ll.il 

Of fundamental importance to our analysis is the conversion of the second-order 
Poisson equation ((T]i) to a coupled system of first-order equations. Introducing the 
velocity vector u = — Vp, equation ([TJ can be rewritten as 



u + Vp = 


in 


(l(t), 


(3a) 


divu = fiS 


in 


Sl(i), 


(3b) 


p = H 


on 


r(t), 


(3c) 


V(r(i)) = u-n 


on 


r(t), 


(3d) 


(1(0) = SI. 






(3e) 



1.2. Some prior results. In the case that fluid is not being injected into the 
Hele-Shaw cell and p = 0, Constantin & Pugh [2] established the stability and ex- 
ponential decay of solutions of ([1]) using the methods of complex analysis. Friedman 
and Reitich [5] also establish this stability result. In [TJ, Chen studied a two-phase 
Hele-Shaw problem with surface tension, and established well-posedness using the 
energy method coupled with certain pointwise estimates from the theory of har- 
monic functions; moreover, he proved that solutions exist for all time if the initial 
interface is a sufficiently small perturbation of equilibrium. Weak solutions have 
been obtained by Elliott & Ockcndon [3] and Gustafsson [7], and classical short- 
time solutions to related problems have been obtained by Escher and Simonett [3] 
in multiple space dimensions. In two dimensions, Escher and Simonett [5] establish 
global existence and stability near spherical shapes using center manifold theory. 

In the case of fluid injection where p > 0, Prokert [8] and Vondcnhoff [9] establish 
global existence results in the case that the injection rate p is a positive constant. 
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1.3. Statement of main results. We propose a simple methodology (equally ap- 
plicable in three space dimensions) for establishing global existence and decay to 
equilibrium for time-dependent injection rates. Unlike the case of zero fluid in- 
jection, the decay to equilibrium is not exponentially fast but, rather, algebraic. 
We use an Arbitrary Eulcrian Lagrangian (ALE) formulation to transform the 
free-boundary problem ((3J to a system of PDE on a fixed domain. This ALE trans- 
formation depends on the signed height function h measuring the signed distance 
between the moving surface T(t) and the expanding sphere B p ^ t y The idea behind 
the proof is the construction of a (total) norm, denoted by ||| • ||| T , which consists of 
a norm || • \\x associated to energy estimates in sufficiently high- regularity Sobolev 
spaces, and a norm || • ||y associated to decay estimates in weaker topologies. The 
total norm is given by 

|||-||| T = sup {\\-\\x+V{t)\\-\\ Y ) 
te[o.T] 

for some function T>(t) — > oo as t — > oo. We shall prove that h satisfies 

|||h||| T < Ce+ \\\h\f T p (4) 

for some integer p. The inequality (j4|) implies that if e ( an upper bound of the initial 
for a norm of the initial data h) is sufficiently small, then |||h||| T stays small for all 
t € [0, T]; a standard continuation argument shows provides global existence and 
the fact T>(t) in front of the lower-order norm || • ||y gives the decay to equilibrium. 
The main results we establish in this paper are the following two theorems. 

Theorem 1.1 (Stability for slow injection). Let (p,Cl(t)) be the solution to ([I]), h 
denote the signed distance between T(t) = dfl(t) and dB p with p(t) defined by p]). 
and the center of mass of the initial domain f2 be the origin. If the injection rate 
p > is such that the corresponding radius p(t) satisfies 

p (fe) (t)(l + t) k pit) 1 

sup — < oo k = 1, 2, and sup -; — < oo for some a < - , (5) 

t>o P(t) t>o (1 + t) Q " 3 w 

then there exists an e > sufficiently small, such that the solution to |T]) exists for 
all time provided that 1 1 ho 1 1 jt 6 (s 1 ) < £• Moreover, the signed distance h decays to 
zero, and 

\\h(t)\\ H -, Hsl) <Cp(ty 2 e-^ Vt>0 (6) 

6ds 



for some constant C > and j3 G (0, 7/8), where d(t) 



Theorem 1.2 (Stability for fast injection). Let (p,£l(t)) be the solution to (JTJ) ^ h 
denote the signed distance between T(t) = dQ(t) and dB p with p(t) defined by ([2]), 
and (xo,yo) be the center of mass of the initial domain fi. Suppose that p > is 
such that the corresponding p(t) satisfies 

/(*)(! + *) . J f_ „ (l + *) a „ 1 3 



sup — < oo and v = sup \ a 



t>o p(f) I t >o p(t) 

as well as one of the following conditions: 

(1) p" < or (2) logp' has small enough total variation. 



sup — < oo > > - (7) 
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Then there exists an e > sufficiently small, such that the solution to (JTJ) exists fi 

r [64v — 21 

all time provided that ||ho|| jjk(§i) < e with K defined by K = max < 6, 



16^-6 



1 



|. Moreover, the signed distance h decays to zero, and 

Vi > 



< c 



VT+t 



Pit) 2 



(8) 



/or some constant C > 0. 

In the following discussion, we define 

f p(t) 2 e^ if 

- 1 p(tf 

I VT+t 

Then ^ and © can be summarized as 

V(t)\\h(t)\\ H2 , {sl} < C 
for some constant C > 0. 



is satisfied, 
if ((7|) is satisfied. 

Vi > 



(9) 



(10) 



Remark 1.3. In Theorem 11.21 t/ie sufficient condition 1 1 h-o II jj fc (S 1 ) ^ 6 which guar- 
antees the existence of the global-in-time solution also suggests that the center of 
mass cannot be too far from the point of injection. 

Remark 1.4. Assumption ([5]) limits the injection rate of the fluid. In this case, 
the location of the center of mass of the initial domain is crucial for the stability 
result, and Theorem \1.1\ states that stability holds when the center of mass of the 
initial domain and the point of injection are identical. 

On the other hand, when ([7]) is valid, the injection is fast enough so that even if 
the center of mass of the initial domain is different from (but close to) the point of 
injection, the stability result holds. 

Remark 1.5. Suppose that ([7]) is satisfied. By Gronwall's inequality, the condition 

p'(t) < Cp(f)(\ + ty 1 for some constant C > 0, 
implies that v < oo. Then for each < £ < 1, there exists C\ and Ci such that 



(i+*r 
p(t) 



< 



i 



p(t) 



> 



i 

C~ 2 



Vt>0 



so that 



+ «)"-= <p(t)<C 2 {l + t)^ 



Vt > 0. 



(11) 



2— i 

Consequently, p(t) iv- < D(t), where we use the notation v to denote v ± e 
whenever < e <C 1. We will make crucial use of the integrability of certain 
functions of these bounds. 



When /j = 0, p(t) = 1 so iJSJ holds. Theorem 11.11 implies that the distance from 
the moving boundary T(t) to the boundary of the equilibrium state (wherein the 
boundary is the unit circle) decays to zero exponentially, which is result that was 
obtained by PQ,[2], [5], and [S], but with more conditions on the data. 
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| ff e (s i) <Cp{t)- 2 Vt>0. 



When p is a positive constant, we have p(t) = yjl + ^t. For this case, Vonden- 
hoff [5] provides a decay estimate for a rescaled function: 

11-11*6(8!) <Cp~ a Vae (0,1) andi>0, 

but this, in fact, shows that h „ 6 , sl . ~ £(!- Q )/ 2 for large t > and hence the 
perturbation may actually grow in time. 

In contrast, for the case of constant injection, our Theorem 11.21 shows that 

||h(0|| H2 .5 (sl) < Cpity 1 Vt>0 

which is the first result proving the decay rate of the actual height function of the 
perturbation of the sphere. With the rescaled variable, 
,,h 
P 

Moreover, our theorem applies to time-dependent injection rates, and thus gener- 
alizes the results of [8] and [9]. 

1.4. Outline. In Section [21 we derive the evolution equation for the signed height 
function h and define the corresponding ALE map ip, which we use to pull-back 
the equations onto the fixed domain B\. In Section [4j we derive some inequalities, 
fundamental for our subsequent analysis. 
In Section [SJ an estimate of the type 

2?(i)||h(i)||^. 5 ( S1 ) < C[||h ||ip.5 (sl) + UNIt + IINIt] V* G [0,T] (12) 

is derived using the decay property of the linearized problem, where [0,T] is the 
time interval of the existence of the solution. We combine this with energy estimates 
for h in a higher-order Sobolev spaces in Section [6l from which we obtain 

—-^\\h{t)\\ 2 HK+1 . Hsl) dt 2 + sup ||h(f)||#K(si)+ p{t)y/ p'(t)||h(i)||j7K-i(si) 

P\ T ) J ie[0,T] L 

< CjM(||ho||fl«(si)) + C a |||h|||^(|||h|||y) + J|||h||| T (13) 

for some polynomial V. Finally, in Section [71 we show how the decay estimate ([12"]) 
together with our energy estimate ([13)) leads to Theorem 11.11 and 11.21 

2. The Arbitrary Eulerian-Lagrangian (ALE) formulation 

We let S 1 = dB\ denote the boundary of the unit ball, and parametrize S 1 using 
the usual angular variable 9. For each 6 G S 1 , let h(9, t) denote the signed distance 
from T(t) to dB(0, p(t)), where the sign of h is taken positive if for x(9,t) £ F(i), 
|x(0,£)| > p(t) and taken negative if |x(6*,i)| < pit). In other words, T(t) can be 
parametrized by the equation 

x(6,t) = (x(6,t),y(6,t)) = (p(t)+h(9,t))N(6) VtfeS 1 , (14) 

where N(6*) = (cos 9, sin 9) is the outward-pointing unit normal to B\. 

Let T(6>) = (— sin(9, cos(9) be the tangent vector on S 1 . Then the outward- 
pointing unit normal at the point (x(6, t), y(9, t)) is 

n(9, t) = 3^(9, t) Up(t) + h(9, t))N(6) - h g (0, t)T(9)} , (15) 



where 3 h (9,t) = ^(p{t) + h(6,t)) 2 + h e (6,t) 2 . 
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Let <I> : M 2 x (0, T) -> K denote the level-set function denned by 



1 r 



$(x, t) = ± |x| 2 - ( p(t) + h( tan^ 1 



Then $ has the property that Y{t) = {(x,y) £ R 2 \$(x,y,t) = 0}. Since the 
boundary is moving with the fluid velocity, $ satisfies the transport equation 

$ t (x(M),t)+u(x(M),*)- (V$)(x(0,i),i) =0 on S 1 x(0,T) 
which implies that h satisfies 

-(p + h)(p' + h f )+ii(x)- [(p + h)N-h e Tl =0 on S 1 x (0,T) 

or after rearrangement, 

h t (6,t)+p'(t) = v(6,t) -Af(e,t) on S 1 x (0,T), (16) 

where v(9, i) = u(x(8,t),t), and 

h e (6,t) 



p(t) + h(6,t) 



T(0). 



(17) 



In addition, we note that with respect to the height function h, the mean cur- 
vature is given by 



H h (x) 



-Qo + h)hge + (p + h) 2 + 2h 2 e 

[(p + h) 2 +h 2 f 2 
-(p + h)h ee +J 2 h + h 2 



I 1 x (0,T). 



(18) 



2.1. A divergence-free velocity. When ho = so that fl = Bi, the solution to 
the Hclc-Shaw equation (|3|) is given by f2(t) = B(0,p(t)), 

p(a; '' ) = ^)-^r l08 ^)' uM = ^w 

In order to have a divergence-free velocity field, we introduce the new variables 
u = u — u and p = p — p, so that ^ is converted to 



u + Vp = in Q(t), 

divu = in f2(t), 

p = H-p on T(t). 



(19a) 
(19b) 
(19c) 



2.2. The ALE formulation. Let ip(-,t) denote the ALE mapping, taking B\ to 
fl(t), defined as the solution to the elliptic equation 

AV> = in B u (20a) 
ij = x on S 1 , (20b) 

where we recall that x = (p + h)N on S 1 . When the perturbation h is close to 
zero, elliptic estimates and the inverse function theorem show that tjj(t) = ip(-,t) 
is a diffcomorphism, and f2(i) = ip{t){B{), By introducing the ALE variables 
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v = uo^i, q=po^i, q = p o ip 7 and A = \7ip 1 , we find that (fl"9|) can be rewritten 
on the fixed domain as 



k^>3 

a;v v ii 



q = H h (x) - q 
h t + p' =v-Af(M) 
where Hh(x) is given by (|T5|) and M is defined in p7|) 



in Bi, 

in Bi, 

on S 1 , 

on §\ 



(21a) 
(21b) 
(21c) 
(21d) 



2.3. The vector JA T N. Let J = det(V# Since i/>(0,t) = (p(t) + h(0, i))N(0) on 
S 1 , we find that 



JA T N 



2 ^ 2 ,1 

"50"' 50". 

(p + h)N - h e T = (p + hyV(0). 







■ ^,1 V> 2 ,2 ' 








'dip 2 dip 1 - 






-V 1 ,! -^,2_ 










he sin + 


(p + h) cos 0, -he 


COS 


f(p 


+ h) sin0] T 



(22) 



2.4. Linearization about the unperturbed state h = 0. When h = 0, h = 
for all t > 0, in which case, ip(x,t) = p(t)x and A = /? _1 Id. Therefore, we may 
decompose (|21a .b) into a linear term and a nonlinear remainder as 



where 



v + p- 1 Vq = / 1 in B lX (0,T), 
divv = / 2 in x (0,T), 



fl = (p-'Sj Af)q y - and / 2 = (^' - pA^ 



(23a) 
(23b) 



In order to determine the linear operator associated to the boundary condition 
(|2Tb). we multiply both sides of ([2lb) by p~ 2 (p + h)" 1 J^, and find that 



1 

q= — 

p 2 



- h ee + p + h + 



where 



1 = 9 2 {P 



p + h 

J6 -P 4 (p + h) 2 



- (7 + l)q-7q, 



(24) 



p 2 (p + h)[j3 + p 2 (p + h)] 
(p + h) 2 (4p 3 h + 6p 2 h 2 + 4ph 3 + h 4 ) + 3(p + h) 4 h 2 + 3(p + h) 2 h 4 + h$ 



p 2 (p + h)[J 3 i + p 2 (p + h)] 



We remark that 



2h 

7 = — 
P 



for some g = C(h 2 + h 2 ,). Since 
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1 ph 



a Taylor expansion in h/p <C 1 shows that q 
writing q = q + p _1 — p'h. where 



2?rp 



p 1 — p'h; thus by 



psh 
2~7rp 



~2n 





h 


;io g (i+^) 






P 



we may write as 



1 



[h ee + h] + p'h - 7 q + £ on S 1 x (0,T), 



(Hat) 



where £ denotes an error-term given by 

hi 



£ 



p 2 (p + h) 



( 7 + l)q + 7p'h- 



P 



2.5. A homogeneous version of (J23J). We can now define a new velocity field 
which is divergence-free. Let /i = w + p _1 Vr, where r is the zero-average solution 
to the following elliptic equation 

p _1 Ar = div/x — / 2 in £>! 

_i <9r 



<9N 



h-N--^- f f a dx on S 1 



(25a) 
(25b) 



We remark that the solvability of (|2"5|) is guaranteed by the solvability condition 



/ (div/x - / 2 ) da; = / f/ x • N - ±- [ f 2 dx) d9. 
Jbi Jo v 2 nJ Bl J 

Let v = v — w, and qr = q — r. By (|25a ) we find that divwi = / 2 ; thus 
that 



u + p _1 Vg< = 
divw = 



9 



hee 

n 2 



p^ p- 

where G = £ — r. We note that on S 1 , 



in Bix(0,T) 
in Bi x (0, T) 

p'h- 7 q + G on § 1 x(0,T) 



implies 

(26a) 
(26b) 

(26c) 



v ■ N = (v - w) ■ M = h t + p' + he(w u T) — (w N) 



, Mw-T) 1 
h * + P H TT T~ 

p + h ZTT 



p + h 



(27) 



3. The total norm and the basic assumptions used for our decay 

estimates 

We first define the total norm ||| • ||| T , used to establish (|4]), as follows: 



||h||| T = 



#K+1.5(gl) 



!,dt 



+ sup 

te[o,T] 



||Mt)||HK( S l)+ p(*)vV 7 W||h(t)|| HK -l (S l ) + P(t)||h(t)|| ff2 .5 (S 1) 



(28) 
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where K = 6 if (0 is satisfied or K = + 1 if J71) is satisfied, and we 

L 16/i — 6 J 

recall that T>(t) is defined in © by 

! p ^2 e 0d(t) if q is satisfied^ 
P(t) 2 ™ 
Fy ' if © is satisfied. 
s/T+t 

We remark that once the boundedness of |||h||| T is established, ||h|| h 2 - 5 ^ 1 ) decays 
to zero at the rate 2?(t) -1 . Throughout the rest of the paper, we assume that p 
satisfies the condition ([5]) or ([7]); that is, 

p( fe )(t)(l + t) k pit) 1 

sup — < oo k = 1, 2, and sup — < oo for some a < - J5j 

t>o p{t) *>o (1 + t) a 3 

or 

+ i f (l + t) a 1 3 _ 

sup — < oo and v = sup <^ a sup — — < oo > > - . (17]) 

t>o p(t) I t >o p{t) > 8 

Moreover, we make the following basic assumption: for t € [0, T] and for sufficiently 
small positive constants e and a to be made precise later, 

\\h(t)\\ H K {sl) +V(t)\\h(t)\\ H 2. Hsl) <a^:l and ||h i| ff K (sl) < e « 1. (29) 

(We will prove that (]2"!)1) indeed holds whenever the initial data is sufficiently small.) 

4. A PRIORI ESTIMATES 

4.1. Estimates of ip. Under assumption (]29j) . elliptic estimates show that 

||W - /9ld|U~ (Bl) < C||W> - pId|| H i.. (Bl ) < CllhHflapi) (30) 
and for 1 < s < K - 1, 

\\D 2 ^\\ H s (Bl) < C||h|| H . + i.. (S i ) . (31) 
Estimate (]30| implies that 

||pA-Id|Uoo (Bl) = ||( P Id-V^)A|| L » (Bl) < C||A|| i « (Bl) ||h|| iI3(s i ) 
< Cp- 1 ||pA-Id|| L » (Bl) + 1 ||h|| 
thus under assumption (]29]) . 

||pA-Id|| L oo (Bl) < Cp^Uhll^^). (32) 
Note that ([30]) and <(32]) together imply that for < s < K - 2, 

\\Dk\\ H , [Bl) <Cp- 2 ||h|| Hs+1 .5 (S1) . (33) 
Furthermore, with J = det(V^), inequality (|30|) implies that 

||J-/» a |Uco (fll) <C||h|| 1 p. (s i ) (34) 
and (]3"2"]) and (]3"3"|) together with (]3"4"]) shows that 

||-DJ|| B »( Bl ) < (7/j||h||jj.+i.5 (s i) , (35) 
from which it follows that for < s < K — 2, 

||D(JA)|| B . (Bl) < C||h|| B , +1 .5 (sl) . (36) 
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log(l + x) - x\ < x 2 if \x\ < 0.5, and 

<9q = 



q = - JL 






h 


log (l - 


--) 




H 2tt 




P> 


9- 



Since 



p'hdh 
p + h ' 



we find that with the basic assumption (j2"9")l . forl<s<K — 2.5, 

||q||H=(si) < CpV 1 ||h||ff ! '(si)||h|| H .( S1 ). 
Since q = q + p -1 — p'h, 



M < C 



P X (! + Uhll/i 2 ^ 1 )!!* 1 !!//^ 1 )) +p'l|h|| ff3 ( S i) 



4.3. Estimates of /i and fa. By the definition of fi and /2, for s > 1, 



<Cp 2 ||h||jj-«+i(si)||q||fl-i.B(B 1 ) + ||h||u-2( S i)||q||j ? ,+i(_ Bl ) 



and 



11/211^-1(230 < Cp 1 ||h|| ff .+i( S i)||v|| H i.B (Bl ) 



\w^)\\v\\h>>(b 1 ) 



(37) 
(38) 

(39) 
(40) 



4.4. Estimates of q. Before proceeding to the estimate of q, we remark that since 

1 



27 -t-s - (7 + l)q + 7P'h - -g , 
p 2 (p + h) p 



hllff^+o-^s 1 ) + P 2 ||h||if s + 2 - 5 (s 1 ) 



in Six(0,T), 



£ = 

by (j3"7| we obtain that 

\\£ ||_h---+°- 5 (s 1 ) < C — ||h|| g 2(gi) 

Note that q satisfies 

Aq = p[div/i - / 2 ] 

q= _h^_h + / h _ 7q + < p Qn §i x (o,t). 
p^ p- 2 

Elliptic estimates together with (|38l) - (f4"Tj) then show that 

!|q!|ff2( Bl ) < C p(||/l||ffl(B 1 )+ \\h\\mB 1 )) + P _2 ||hj| ff 3.5 (S l ) + p'||h|| H 1.5(Sl) 

+ Il7q||ff i - 5 (s i )+ ll^llffi-^s 1 ) 

< C ||h|| ff 2 (S i ) (||q|| fl -2 (Sl) + ||v|| H i (i3l) )+p- 2 ||h|| / j3.5 ( si ) +p'||h|| Jf i.5 ( s 
+ l|h||/f 2 (s 1 )l|h||ff 1 - 5 (s 1 )+ P _2 H n llH 2 (s 1 )ll n llH 2 - 5 (s 1 )+ ll7q||i? 1 - B (S 1 ) 



(41) 

(42a) 
(42b) 



Since 



||v||tfi (Bl ) < p 1 ||q|| jF f2( J B 1 ) + 11/illffMBi) < C (P 1 + ll h llff 2 (S 1 ))l|q||if 2 (Bi) ! 



l7q||/fi-=(si) < c 



p- X ||h| 



P" 2 ||h| 



q||/j 2 (Bi), 



by assumption (f2T)|) we find that 



qllff 2 ^!) < C p 2 ||h|| J? 3.B( S l) +p'||h|| H 1.5( S l) 
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Similarly, elliptic estimates also show that for integers 1 < k < K — 1, 

||q||^+i( Bl ) < C ^ 1 ||h|| H&+ i ( si ) (||q|| H i.5 (Bl) + ||v|| ff i.5 (Bl) ) 

+ ^ 1 |ih|| ff 2 (S i)(||q|| ff A ; +i( Sl ) + ||v|| H fc( Bl )) 

+ /5 _2 ||h||_H-fc+2-5(si) + p'||h|| fl -*+o.5( S i) + ||h|| ff 2(si)||h|| i jfe( S i) 

+ P~ 2 |l h lli? 2 (s 1 )ll h llff )e + 1 - 5 (S 1 ) + C||7q|U fc +°- 5 (s 1 )- 



Since 



|v|| H a (Bl) < p iqlUa+i^) + ||/i||h=( Bi ) 

< Cp~ l \\(\\\ H s+i [Bl] + ||h||H3.B(si)||h||H«+i( S i 



and 

||7q|| ffs +o. 5( gi ) < C 
< C 

we conclude that 



ff 1 (§ 1 )l|q|l# s + 1 (-Bi) + ||7l|iT«+ - B (S 1 )ll < ll|.ff 1 - fi (B 1 ) 

fp(§i))l|q||.ff»+ 1 (-Bi) 



(^HhH^^+^llh" 2 



(43) 



+ P INI-H" S +0.5(§1) + P 2 ||h|| /i -2(gl)||h|| ffs + 1.5(Sl)||h|| iJ -3.5(sl) 



l|q||ff"+i(Bi) ^ c [p _2 |Nlff*+ 1 (S 1 )Nlir !J (B I ) + P~ 1 |l h llff 3 - 5 (s 1 )ll h llff fc + 1 (s 1 ) 

+ /5 _2 ||h||_H-fc+2-5( S i) + p'||h|| JT *+o.5( S i) + ||h|| ff 2(si)||h|| i jfe( S i) 

+ P _2 ||h||ff 2 (s 1 )ll h lli/ fc + 1 - 5 (s 1 ) 

In particular, when k — 2, using our basic assumption (|29[) . we obtain that 

||q||ff3 (Bl ) <C p^WhWH^iS^WH^^) 

+ P~ 2 \\Hh*h§i) + p'l|h||i/2.5 ( si) + llhU^^i) 

+ P _2 ||h||ff2(sl)||h|| ff 3.5(Sl) 

and is hence made small by (|29j) . This, in turn, implies that 
ll v llff fc (Bi) + P~ 1 \\ c l\\H k + 1 (B 1 ) < C p~ 2 ||h|| ff 4.5(si)||h|| fl -&+i ( gi ) + p~ 3 \\h\\ H k+2. 5{s 



+ — l|h|| ff fc+o.5(si) + p ||h|| ff 2( S i)||h|| H &+i.5(si) 

< Cp^ 1 ||h|| i? fc + 2.5( S l) + p'||h|| JJ *+0.6(Sl) 

Consequently, 

II/i||h s (Bi) + ||/2||^-i(Si) 



(44) 



< Cp~ 



||h||ff»+2.5(Sl) + p'UhH^f^O.SfSl) 



(45) 



for 1 < s < K — 2, where we have used the basic assumption ||h||#-2.5(§i) < a. 
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4.5. Estimates of r and w. Applying elliptic estimates to (|25|). using (|45[) we find 
that 



||r||ffs + i( Bl) < Cp 1 ll/illfl-^so + ll/allfl-'-U-Bi) 

< Cp' 3 HhHiya.s^i) ||h|| ffs + 3 . 5( si) +p'||h|| H a + o.5 (s i ) 
Since w = /i — Vr, we also obtain that 

||wj| ffs(Bl) < ||/i|| B =( Bl ) + ||r||ff3+i( Bl ) 

< Cp~ ||h|| B 2.5(gl) ||h|| ff .+2.5( S l) + p / ||h|| H -a+0.5(§l) 

4.6. Estimates of v and g in terms of h. Since v = v — w and = q — r, 

< l|v|| ffs(Bl ) + ||w|| ffs(Bl) +p~ 1 ||q|| H3 +i(s 1 ) +P _1 ||r!lff=+i(B 1 ) 



(46) 



(47) 



(48) 



(49) 



< Cp 1 ||h|| i? «+2.5(gl) +p'||h|| H s+0.5(Sl) . 

Moreover, since Vg = pv = pv — pw, 

IMI.fp+i(.b 1 ) < p[||v|| B s( Bl ) + ||w||ir.( Bl )] + \\q\\m(B 1 ) 

< Pll V llH-(Bl) + C |jh|| ff 2.5 (S l) + / 9'||h|| J J0.5( S l) 

which together with (|46|) further implies that 

llqll^ + ^Bi) < P\W\\HHBt) + C ||h|| H 2.5( S i) + p'||h|| B o.5 (s i) 

+ Cp _3 ||h|| H 2.5(gl ) ||h|| H i+2.5(Sl) + p'||h|| ff3 + 0.5(Sl) . 

Therefore, 
||q||iy-+i(Si) 

< /f||v||i? S (B 1 ) + C(l + p') ||h|| B 2.5(Sl) + p _3 ||h|| ff 2.5 (S l)||h|| B3 + 2.5(Sl) 

On the other hand, we also have 

l|Vq|| ff . (Bl ) < ||(pA T -Id)Vq|| 

< p\\v\\h'{b{)+ C ||Id - pA|| L oo ( si ) ||Vq|| ffs(Bl) + /9||DA|| WJ .- I ,4 (Bl) ||Vq||jr4( Bl ) ; 

thus by estimates (|32|) . ([33]) and assumption ([29]) we find that for 1 < s < K — 1, 

l|Vq|| B .,+i (Bl) <p\\Ah-(b 1 )- (50) 

4.7. Elliptic estimates for the height function h. We may rewrite the bound- 
ary condition (|26b ) as 

-h.ee + P 2 p'h = h + p 2 [q + 7 q - S] . 
Elliptic estimates then imply that for 1 < s < K — 1, 

||h||!r.+2.B (s i ) +P 2 p'l|h|| B3+1 .5 ( si) 

< c||h||^ + o.5 ( si) + p 4 [l|q|| BS +i( Bl) + ll7q|| BS+ i (Bl ) + 11^ IIh-+o.«(s») 
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thus by (|4Tj) . (|43|) . (|46|) and (|49| together with the smallness of ||h||#e(§i), as well 
as the boundedness of p'/p, we find that 

\\h\\ 2 Hs+2 . Hsl) + p 2 p'||h||^ + , 5(sl) < C[p 6 ||v||^ (si) + ||h||^ + , 5(sl) 

+ P 2 ||h||H S+ o.5 (S i) + P 4 (l + P' 2 )l|h||^.5 ( si) . 
It then follows that for 1 < s < K — 1. 



P I|h|| ff3 + 2.5(S1) < C 



.-3 



(51) 



+ PP' 2 )l|h||^ 2 . 5(s i ) 

4.8. An elliptic estimate via the Hodge decomposition. 

Lemma 4.1. Suppose that v G L 2 (B\) satisfies divw = curlw = 0, and the tan- 
gential derivatives d e v € L 2 (Bi) for all £ = 1, 2 • • • , k. Then v € H k (Bi) and 
satisfies 



\h»{b 1 ) < \\9 3 v\\ L 2 {Bl) 



(52) 



Proof. For k e N, 

IMU*(.Bi) < C \\v\\ L 2 (Bl) + ||curlt;|| i?i ;-i (Bl )+ ||divw|| H(i -i (Bl) + \\v ■ N|| J? fe+o.5 ( si) 
Since curlu = divu = 0, we find that 
IMIfl*(Bi) <C \\v\\ L 2 (Bl) + \\v ■ N|| fffc -o.5( S i) 

< C 

< C „ 

L 2 (S!)+ \\d 2 V ■ N|| if k-2.5( S 1) + ||<% • T||_ ff fe-2.5( S l)+ ■ T)||^*-2.B( S l) 

l 2 (Bi)+ ll<9« ' N lli/ fc - 2 -5(si)+ ||9 2 w • N|| H fc-».5 (s i)+ \\Bv ■ T|| H -fc-3.5 (s i) 



< c 

< C 



< C 



l 2 (Si)+ l|5(w • N)|| fffc -i.5(si) 

L 2 (Bl) + \\8V ■ N|| iffc -1.5 (S l ) + \\V ■ T|| fffc -1.5 ( Sl) 



L 2 (Bi)+ E ' N llff- 5 (S 1 )+ W dv ' T llff-°- 5 (S 1 ) 

3=1 



Again, with divw = curry = 0, we see that ||u ■ N|| 5 (§i) < C|| t'll fc (Si ) 1— ' 

5. Decay estimates for ||h|| ff 2.5( S i) 

5.1. The Fourier representation of solutions. Let Gi = G + 7q, and express 
h and Gi in terms of their Fourier series: 



h ^ *) = 4=" E M*)^. Gi((9, *) = i E Gl * (*) 



Since g is harmonic in i?i with the Dirichlct boundary condition (|26b ) , the Poisson 
integral formula shows that 

^2tt 



?(r, 0, t) = J= V [^^h fc + p'h fe + G u 



for r < 1 . 
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Taking the inner product of (|26|) and the vector Af, by ((27)) and that 



we find that 



dq dq 



h t + p' + -^- 



h e Vq • T 



p dr r=i p{p + h) 



h g (w ■ T) 



1 



\ \k\m 2 -D . y 



E 



r=l p + h 

h fe e^ = - 



+ 



dq 



h s (w ■ T) 



p(p + h) 89 
Therefore, for k ^ 0,±1 we have 



r=l 



p + h 



2irp 



■=i + 2tt 



(53) 



/ 2 da;. 



dhfe 



l*l(l*l 2 -l) 



h fc = 



(54) 



where Rf. is the Fourier coefficient of R defined as the right-hand side of (f5"3")l . Note 
that estimates (|4T))) and (|46|) then suggest that 

ll-Rllff^^S 1 ) < C ||G + 7q|| ff 3.5 (s i) +p _1 ||h|| ff 2.5 ( si ) ||q|| ff 4 (Bl) 

+ P _1 ||h||/i-3-5(gl)||q|| ff 2.5( Bl ) 

< C i o~ 1 |jh|| ff 2.5 ( si ) p'^lhllija.s^i) + / 9 / ||h|| H 3.5( S i) . (55) 



5.2. The if 2 5 -decay estimate. Define 



«t) 



: l*l'-') t | t |4y 



|fc|(|fc| 2 -l) ^ 



P(i') 3 



df + |A| logp(t) . 



Then the use of I k (i) as the integrating factor in (|54"]) implies that 



M*)l < 



Since d(t) + 2\og p(t) = I 2 (t) < I k {t) for all \k\ > 2, and p(s) < p(t) for all s < t, 
we find that 

|h 0fc |+ / p(s) 2 e d ^\R k (s)\ds V|fc|>2; 



p(t)= _ 

thus by Holder's inequality, 

J2 (i + ifci 5 )ih fc (t)i 2 



fe#0,±l 



< 



Ce -2d(t) 



p(ty 



l|ho||^. 5(SI) +i / p(s)*e 2d M\\R(s)\\ 2 H2 . 5{§1) ds 



(56) 



Since p'(t) is bounded, by (1551) and interpolation we obtain that 



p( S ) 4 e 2rf ( s )|| J R( S )|| 2 f2 . 5(sl) d S 



(57) 
ds. 
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5.2.1. The case that p satisfies ((5]). In this case, by interpolation ()57[) implies that 

+ c^( s ) 2 p'( s ) 2 ||h(.s)||^ (sl) !|h( s )||^ 5 5 (sl) d s . 

Since ||h.(t)|| ff e (s i ) + p(t) 2 e^ ( * ) ||h(t)|| H 2.5( S i ) < |||h||| T , we find that 
p( S ) 4 e^\\R( S )\\l 2 ,, {sl) ds 



< C 



< C 



+ 



V p( S ) 32 /7 ^ p( S )34/7 
p'( S f 



ds 



Jo Vp( S )"/7 - p (s) 13/7 

thus by ([5]), we obtain that 



+ pwy2-? ms ) d , {s) \ ds 



\\T' 



p(s) 4 e 2d ^\\R(s)\\l 2 . Hsl) ds<C / eP-WMd'Wds 

Jo 

C 



o (2-f(3)d(s) 



2-fP 



s=0 



16 7 

if 2 /3 > or /3 < - . Inequality (j51)f then implies that 

7 8 

„-2d(t) r 

£ (1 + |*| 5 )|h fe (t)| 2 < C^j- fllholl^-op) + CflteC-^^ffllill* 



\k\>2 



< c 



e ~2pd(t) f 



||h |||, a .» ( si ) +^|||h|||; 



(58) 



t>o 



where we conclude the last inequality from supte 2 ^*- 1 / 7 < oo which is a direct 
consequence of assumption ([5]) as well. 

5.2.2. The case that p satisfies assumption (|7|). In this case, 



d(t) 



6 



Jo Pi*) 3 
thus d is bounded. Since 



'o (1 + *) 3 "" 



-ds < oo V t > 0; 



||h(*)llff»(Si) +P(tf ^\Mt)\\H*H&) < IIN| T , 
by interpolation, (|57|) implies that 



2+2 K ~ 5 - 5 

I 1 K — 2 . 5 



+ C I p( S ) 2 p'( S ) 2 ||h( S )|| 2 / 2 2 ™||h( S )||™ )( i S 



< c 



l-4u~ 2K-8 



1-Av~ 2K-6 



p(s) ^ K ~ 2 - 5 + p(s) K ~ 2 - 5+ p'(s) 2 )ds 



16 



Note that by the definition of K, K > 
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64z/ - 21 



16v - 6 



; thus by ©, 



l-4u~ 2K-8 



p { s )^7^^2T, + p( s ) i ^ra+ 2 (0 '( s ) 2 < C a {\ + t)- a 



for some a = a(v) > 1. Therefore 

te -2d{t) ft 



R(s)\\H* S(§i) ds < 



C a t 



{\ + t)- a ds < 



Ct 
W 4 



and we then conclude from ([55| that 
Pit) 



f-L- Y, (1 + |fc| 5 )|h fe Wl 2 < C[|!h ||^, (sl) + |||h||£ 

|fc|>2 



(59) 



5.3. The mass and the center of mass. For x G we write x = (x,y). 

Equation ([2|) shows that J^r t \ — ^p(t) 2 ; hence, 



irp(tf 



(xdy — ydx) 



2jv 



( P (t) + h(6,t)) 2 d8. 



As a consequence, 

h (t) = -ih\(i) = -^£h fe h_ fe 



2p 
1 

2^ 



Ihnl 2 



1 fc/-l,0,l 



(60) 



Moreover, by Green's identity and the fact that that / HndS = 0, we see that 

'r(t) 



(x,y)cM 



n(t) 



r(i 



(x,y){u ■ n)dS = - I {x,y)—dS 



r(t) 



d(x,y) 
dn 



pdS 



dp 

r(t) 
(x, y)Apdx 



n(t) 



HndS* = 0. 



r(t) 



Therefore, the center of mass of f2(i) does not change in time; thus 



and 



xo 



yo 



/ xdA= / r 2 cos9drd9 = - {p + hfcosOdO 
JQ(t) Jo Jo 3 Jo 



ycL4 = 



fi(t) 



27T /.p+h 



JO 



1 

3 7o 



Letting xo + iyo = r o e > we see that 

r2w 



sin (9drd(9 = - / (p + h) 3 sin 6W0. 



(p(t) + h{9, t)fe ±l9 d6 = 3r e ±ie °; 



o 



thus 



p 2 h ±1 +ph 2 ±1 + -h3 ±1 =3r e ±l 
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In particular, since fgj = fkgj-k, we find that 



hi = — 2_, h fc hi_ fc - — ^ h fc _<h < h 1 . 



3r e ie ° r2h , l|h||£ 2(sl) 2|h | 



L o 



3p 2 



3p 2 



h x -(-+ — ^ — )E h ^ 



H \t~2\>2 H fe^0,l,2 \l-k\>2 



(61) 



By (j58[l and the Schwartz inequality, 



^ h fc h ± i_ fe < 

fc^O.zbl 



< 



D(t) 
C 



|ho||//2. 5 (gi) + |||h||| T 



V jhfcj 
^ 1 + Ifcl 2 - 5 



|ho| 



r]lNU»(S 



(62) 



On the other hand, since \hj\ < ||h||i2(§i), 



^ h 2 -£h f h_! < 

\e-2\>2 



< 



C 

W) 
c 

W) 

and similarly, 

fc#0,l,2 |£-fc|>2 

C* 



|ho||^ (S1) + |||h|||^ £ ~ 2 |L 6 



\l-2\>2 



|h || ff2 .5 (Sl) + |||h|||^||h| 



2 

L 2 (§!)l 



(63) 



< 



< 



v{ty 
c 



V{tf 



h 0|lH2.5 (S l) + |||h|||yj ^2 ^2 

fc#0,l,2 \£-k\>2 



(l+|fc-£| 2 - 5 )(l + |fc| 2 - 5 ) 



h o||^-5(Si 



|L 2 (S! 



(64) 



Moreover, by assumption (|2^|) , |h | < 2Tr<r <C 1 and ||h||x,2(si) < \Fhi(J < 1. As a 
consequence, (f6"Tj) together with (|62|) . (|63| and (|64"]) implies that 



|hi|< 



C 



V{t) 



|ho||if2-5(si)+ |||h||| 2 -,+ |||h|||y 



(65) 



A similar argument also suggests that |h_i| shares the same upper bound as |hi|. 
Therefore, once again using the inequality |h | < 2-kg <C 1, (|60]l together with (|62|) 
implies that 

|ho| < Cp- 1 [|£i£_i| + 7^y(l|ho|| ff 2 . 5(sl) + llhlll^llhll^sx)" 



< 



C 
©(*) 



|ho||ff2.5(§i)- 



||h||| T + |||h||| T 



(66) 
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Combining flM} (or ([59])). ([65]) and (|66l), we conclude that 



[£(l + |£f)|h 



< 



C 



2>(t) 

p(t)- 3 + p(t)p^ft) 2 



0|l_ff2.5(gl)- 



dt < oo 



(67) 



Furthermore, since 



for 2? in both the slow and fast injection cases, (jSlj) implies that 
p(t)- 3 ||h(i)||^ K+1 . 5(sl) * 



< C 



p{tf\\v(t)f H ^ 1{Bi) dt + C ||h ||^.5 (sl) + |||h|||* + 



16 



(68) 



6. Energy estimates in the space of higher regularity 
d 

Let d = T- V = — denote the tangential derivative. Tangentially differentiating 
88 

pTk ) f-times (£ = 0,1,2,3,4,5) and then testing the resulting equation against 
Jd e v, we find that 



J\d £ v\ 2 dx 



JA^q^dVdx 



E 

fc=0 



Bi 



Jd e - k Aid k ci, j d e v i dx. 



Writing the second integral on the left-hand side as 



JAid^j dVdz = / JA?'(^q),i<9Vd:r + / JA^ [^q, J -(^q), J ]aVda; 



Bi 



then integrating by parts in Xj leads to 



/ J|9 £ v| 2 dx+ / JA^qSVN^S = - V ( f ) / 

is! js 1 fc=0 V™/ 



JA? 



Bi 



(69) 



Using the Kronccker delta symbol Soi, which vanishes for all I ^ 0, by (|33|) . ([34]) 
and ([5Cl|) together with the continuous embedding if {B\) L p '(Bi) we find that 



JA? 



[B £ qiJ -(9 £ q)„ ] dV + <9 £ q[(<9V)„ -9 £ v 



dx 



< Cp 2 (l - 6 e)\\v\\ Hl -i {Bl) \\v\\ He{Bl) , 
as well as for < k < £ — 2, 



Bi 



Jd*- k Ald k ci, j d i v'>dx <<V||DA|| 



ollVqll 



W k - i (B 1 



\H'(B 1 



< Cpllhlii/f-fc + ifsijIlvllfl-fc+o.s^^Hvljfl-^B^ 

and for k = £ — 1 , 

j^- fc A^ fe q„ aVdxj < Cp 3 ||DA|| L » (Bl) ||v||^-i (Bl) ||v||^ (Bl) 
< Cpll ^s^gi ) || v|| (Si) H v|| ^ t . 
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Therefore, (|44|) suggests that 

J\8 e v\ 2 dx+ [ JA^qflVNj-dS 



S 1 



< / JAi^q^V !3 -da: + C(l-^)p 2 ||v||^- 1(Bl) ||v||^ (Bl) 
+ C(l - <5o«)||h|| ff £+i( S i)!|vj| ff f( Bl )||v|| fl -i( Bl ). 



(70) 



e-i 



By flUb), A;//v'. ; = ( l )d l - k Aid k v\ f , thus for t > 1, 

J^qA^Vyete < Cp 2 (l - M[l|Vq||^-i( Bl) P A IU»(Bi)ll v ll^(Bi) 



Bi 



<-2 



+ ||Vq||.ff/-o.5( Bl ) E ll- DA llB"^ fc - - 5 (Bi)ll v llH fc + 1 (-Bi) 



fc=0 



< Cp(l-S oe ) ||h||^+i(si)||v|| H i( Bl )||v|| ff i( Bl )+ ||h|| ff 3(gi ) ||v|| ff £-i (Bl) ||v||^ (Bl) 



We now focus on the second integral of the left-hand side of (|70|) . By identity 
22|) and the boundary condition ()21b). 



JA^qaVfydS = / ( P + h)(d e vAf)d e qdS 



i-x 



s 1 



(p + h) a £ (v • AA) - ^ 9 fc v • B e qdS 



k=0 



[ ( P + h)aV • N)B e \ - j h 3 [(p + h)h ge ] + p'h 



+ / (p + h)5 e (vM)d e J h 3 (4 + h2)-p- 1 +q 
Js 1 l 
t-\ 

E 

fc=0 



<2S 



(p + h)<9 fc v • d^Md^qdS. 



By employing the i?°- 5 (S 1 )-7J 5 (S 1 ) duality pairing, 



53 / (p + h^v-^'WaVS <C(l-^)p[||v|| L 4 (sl) ||^AA|| i2(sl) ||^q|| L 4 (s 



fc=0 



+ (1 - (5 f - <5w)||<9^ 1 v||i/-o.5 ( si ) ||^q||/ / o.5 (S i ) ||5A r || ff i(si) 



+ (1 - 5 oe - 8 U - 5 2 e) l|S fc v|| H -o.5 (S i ) ||a'q|| H o.« ( si ) ||a < k Af\\ H H§i) 
fc=i 

< Cp 2 \\v\\hHb 1 )\W\\ h i(b 1 ) [l + P _1 ||h||jj«+i(si)] 

+ llvll^-ifB^HvllH^BO [l + /0 _1 ||h||H4 (s i)] 
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and by dSTJ), 



(p + h)d e (v ■ j\f)d e cidS < Cp\\v ■ Af|| ff £-o.5 (S i)||q|| H <+o.5 (S i) 
< Cp || v • A/"|| fl -f-o.5(si)||h|| ff 2(gi)||h|| ff t+o.5( S i). 



Moreover, since 



2ph + h 2 + h 2 , 



and 



we find that 



J h p(p + Jh) 

2(p + h) 1 rh(2p + h) + h 2 h 2 + h 2 



Jh(p + Jh) p L Jh(p + Jh) 2 Jhp(p + Jh) 



( P + h)d e (vM)B )! [J^ 1 -p- 1 }dS= / (/? + h)^(v-AA)9 



h 2 + h 2 



< 



+ / ( P + h)d t (-v-M)9 
2(p + h) 



s 1 
dS 



2h 



Jh(p + Jh) 



dS 



Jh(p + Jh) 



J h p(/9 + J h ) 

^(v-AA)a £ h^ + Cp- 2 ||h|| H2 .5 (sl) ||h||^ +1 .5 (sl) ||v||^ (Bl) 



1 



<- / d*(v -Af^hdS + Cp 2 ||h|| i? 2. 5(s i ) ||h||^ + i.5 ( si ) ||v-A/'||^-o. 5( si ) 
P Js 1 

+ Cp _2 ||h||H2.5 ( si ) ||h|| iJ -«+i.5(gi ) ||v|| ff f (Bl ). 
Finally, because of the identity 

(p + h) 2 _ 1 (p + h) 4 (2p + h)h + 3(p + h) 4 h 2 + 3(p + h) 2 h 4 g + h« 

J h p~ j h/3 [j h + P (p + h) 2 ] 

by the evolution equation ([2TH) we obtain that 

(p + h)0*(v • AQ^ [ - J h 3 (p + h)h ee + p'h] dS 

J^ 3 (p + h) 2 d e {vAf)d e+2 hdS+ f (p + h)p'8 e {v -Af)d e hdS 

Js 1 

E (I) X + h )^( v • ^~ fc [ J h 3 (^ + *)] ^ fc+1 h^ 

(p + h)p'a f (v-AA)a £ hdS* 

- C||h||/J-2( S 1)||V • Af\\ H e-0.5( S l) p _2 ||h|| ff 4 + 2.5( S l) + p'||h|| fl -* + 0.5(Sl) 

- Cp" 1 ^ ■ ■A/ r ||i?«- ( '- 5 (s 1 ) l|h||//2(si)||h||^ + o.5( S i) + ||h||_ H - 3 ( S i)||h||. ff i+i.5( S i) 
Since div(JA T v) = 0, by the normal trace estimate, we find that 

\\vAMh-°-*<v) = |Kp + h)- 1 JA T v-N|| H -o.5 ( s I ) 

< (p + h)- 1 JA T v|| ia(fll) + ||V(p + h)- 1 JA T v|U2 (i3l) ] < C||v|| L 2 (Bl) 
and similarly, for < £ < K — 1, 

||v-JV|| jr t-o.. (a i) < C(l + p- 1 |||h||| T )||v|| H . (Bl) . 



k=0 

> — — \\A e + 1 h\\ 2 
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Therefore, (|70]) implies that 



jf P J\d*v\ 2 dx + IA||^ +1 h||f 2(sl) + jf pp'( P + h)^(v • AO&hdS 



< / a £ (vAA)^h^ + C(l-^)p 3 ||v||^- 1(Bl) ||v||^ (Sl) 



c 



P 1 ||h||_H-«+2.5 ( gi) + pp'\\h\\ H 1+0.5 (^i) Ijhlj^.Sfs^Hvlj^f^). 



By Young's inequality, the inequality above further implies that 
d 



p d Wv\\i HBl) + ^ll^ +1 h||i 2(s i) + 2 J pp\p + \x)<f{vN)^dS 

<2 f d e ^-M)d e hdS + C(l-S 0( )p^v\\l t - 1(Bi} +Sp^v\\ 2 Ht{Bi} (71) 
+ Csp 1 p' 2 ||h||^-2.5(si)||h||^ +1 (g 1 ^ + C$p J ||h|j^-2.5(si)||h||^ + 2.5( S i)- 
6.1. The case £ = 0. Using the iJ - 5 (§ 1 )-ff-°- 5 (§ 1 ) duality pairing, 

/ ( PP '(p + h)(yAf)hdS + / (y-/J)hdS <C(p 2 p' + l)||v|| i2 ( Bl )||h|| H o. 5(s i) 
< ^(pp' 2 + p- 3 )||h||^ 2 . 5(sl) + 5p 3 ||v||2 2(jB0 . 

1 p'(s) 2 



Since 



ds < oo , 



lp{s) 3 V(s) 2 p(s)V(s) 4 . 
choosing S > small enough and integrating in time of (|71|) over the time interval 
(0,t), by §7$) we find that 

ll<9h(t)|| 2 L2(s i) + f p( S ) 2 \\v( S )\\ 2 L2{Bi} ds 



10 



<C / p( S )p'(s)ih( s )||^. 5(sl) d s + C ||h ||^. S(S1) +|||h|||*+|||h||£ 



Define 



7V(||ho|| H K (S i), |||h||| T ) = ||ho|| 2 ff K (s i)+ |||h|||*7>(|||h||£ ) (72) 
for some polynomial function V . If ([5]) is satisfied, due to the exponential decay of 



#2.5(§), it is easy to see that 
ft 



p(s)p'(sy\\h(s)\\i I2 . Hsl) ds < C7V(||ho|| ff K( S i), |||h||| T ). 
Now suppose that (J7J) is satisfied. Then ([TT]) and (j6"7| imply that 



p(s)p'(s) 2 ||h(s)||^. 5(s i)ds<C 



p(*) lA 



-ds 



Ar(||ho|| H K (S i),|||h||| T ) 



< C 



(! + *)■ 



/o (l + 2+ - 
thus in either case, 

r' 



o (l + *) 2 p(<0 
iV(||ho||HK( S i), |||h||| T ) < CW(||h || ff K (S i), |||h||| T ); 



l|5h(t)||i 2(S i) + / p( S ) d ||v(s)||i 2(Bi )ds < C7V(||ho|| ff K( S i), |||h||| T ). (73) 
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6.2. The case 1 < I < K - 1. Define h = h H . Then 

2p 

P 2p 2p p 

Therefore, since ||v|| ff i.5 (Sl) < Cp^ 1 ||h|| ff 4 (s i ) + / o'||h|| ff 2 (s i ; 



by (gU), we find that 



pp'(j> + h)5 t (y-Af)d e hdS- - / p'{p + h)8 e {v-Ar)d e (h 2 )dS 



> / pp'd e [(p + h){v-N)]d e hdS - V (?) / pp'd e - k hd k (vAf)d*hdS 



- Cpp'l|h||ff2(si)||h|| H « + i (s i)||v|| ff f( Bl) 

> / pV'^(h t + ^h)^h^-c 5P -V' 2 !|h||^. 5(sl) !|h||^ +1(sl) 



p 

3|U,I|2 



(74) 



- V ||v|| ff£(Bi) - Cp'^||h|| ff 4 (S i) +p / ||h|| ff 2 ( si 
Since Pv / P 7 |l ll || jF/ 5 (s i ) < |||h||| T , 

p'||h(.s)|| ff 4 (sl) + p'( s f\\h( S )\\ HH§1) } |j^h||| 2(sl) d S 



d e h\\l 2{§1 



< 



p'||h(s)|| ff 4 (S l ) +p'(s) 2 ||h(s)|| ff 2 (S1) 





r°° I 
'/ -72 



< c 



t P '||h(.s)|| H 4 (S1) +p'( S ) 2 ||h( S )|| ff2(Sl) ||^(h 2 )||| 2(sl) d S 

o P L 

h(s)|| ff 4 (s i) + p'(s) || h(s) lU^gi) ■ 



P ( S )2 

+ C7V(||ho||HK( S1 ),|||h||| T ) 
By interpolation, 

||h(t)|| ff 4 (Sl) < C||h||™ ^Uhll™) < CP(t)-(™)|||h||| T . 
If ([5]) is satisfied, due to the exponential decay it is easy to see that 

IhWii^ + p'WiihWjU^^^ 



(75) 



p{sf 

Suppose that (J7J is satisfied. Then 

) ||h( S )|| H 4 (S1) + P '(s)||h( S )||^(S 1 



(76) 



f 

Jo 



Pis) 



■ ds 



< C 

< C 



1 



o V(s) 2+(2 ~^ )( ™ ) (l + *)p(«) 

1 



ds 



o (1 + s ) 2 ^+( 2 ^-5)(K i Sir) 



h||| T + C|||h||| 5 
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By the definition of K, the exponent of the integrand 2v + (2v — — ) (— — — — ) > 1 
if v~ > 3/8, so j76]| is still valid if © holds. Therefore, (|75) implies that 



p'||h(s)|| i? 4 (s i ) + p'(s) 2 ||h(s)|| ff 2 (S i 



\d e h\\l 2{sl) ds 



< C\\\h\f T + CiV(||h || ff K (sl) , |||h||| T ) < C 5l iV(||ho||HK (sl)! |||h||| T ) + <*i|||h||£. 



Since 



f^p , d t (h t + —h)d e hdS- 
P 



1 d 
2di 



p 2 p"\\d £ Hh { ^ 



the combination of ((711), (EJ), (JT5J) and (JTHJ) suggests that 



ll^ +1 h|| 2 L2(sl) +p 2 p'll^h| 



U<p)+ P S Wv\\l. {Bl) ds- p'p"\\d e M\i^)ds 
Jo Jo 

< ll^h||| 2(sl) + C Wl iV(||ho!| ff K (sl) , |||h||| T ) 

+ Cs So p3 " v "ff f - 1 ( B i) rfs + a / /° 3 ll v llH«(Bi) d * + *illl h lllr 

for some constant < c < 1. 

On the other hand, for < I < K — 1 , assumption (j2"5)) and condition (JSJ or ([7]) 
imply that 



pV ||^h||£ 2(§1) < 2pV 



|^h|| 2 2(S1) + ^||^(h 



2\l|2 

L 2 (S! 



<2p 2 p'||a f h" 2 



Q^|h||^ 1(S1) ||^h|| 2 2(sl 



< 



2pV||^h||£ 2(§1) +C*Mt 



which, in turn, implies that 



\\B £+1 h\ 



l 



L 2 (S! 



p 2 p'\\d e h\\l 2{sl) + I p*\\d\\\l HBi) ds I p*p"\\d*h\\l Hsl) ds 



„3n m. ,n2 



„2 "i 



< ll^h|| 2 2(S1) + C«, ffl JV(||ho||flK( 8 i), IINIt) 



(77) 



c 5 



p 3 ii v ii^- 1 (si) rfs + 5 / p z \\Am {Bl) ds +^im?T 

<J 



ii Ca is chosen to be smaller than 2S±. 
If I = 1, by CHI) we find that 



|9 2 h||i 2(sl 



Vp'IW 2 



L 2 (S! 



P 3 ||9v|| 



i 2 (S] 



\ds 



p 2 p"||<9h|| 2 L2(S1) d S 



< C^iVdlhollffK^i), |||h||| T ) + 5 p 3 \n\ 2 HHBl) ds + 2<5 1 |||h| 



which, combined with Lemma 14.1 1 and (|73|) . implies that 



pVllh" 2 



c / P \W\\m(B 1 ) ds - P p"\\ h \\m(si) ds 



<C7 il 7V(||h || ff K (sl) ,|||h||| T ) + 2 ( 5 1 |||h||| 2 . 
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for some constant < c < 1. A similar argument, again relying on Lemma 14. 1[ 
shows that 

||h||^ K(s i ) + ip 2 p'l|h|lHK- 1(s i ) + 

< C Sl N(\\h \\ HK{sl) , |||h||| T ) + <75i|||h||| 2 . 

Finally, we look for an upper bound of / ( 2 ( 0"||h|| 2 ^ K -i ( -gi' ) <is to close the en- 



P ll v llffK- 1(Bl) rfs- / pV'HhH^K-ipijds 



i(Si) 1 

*/ U 

ergy estimates. If © is satisfied, by interpolation and the inequality ||h||# 
i 



6( S 1) 



< 



p( s y P '(s)\\h( s )\\ 2 HHS1) d S 
°° P ( S ) 13 / 7 



< c 



1 + s 



< C 



I + s ii-v-/nj 1 —(Si)ll"WllH e (S 1 ) 1 

7V(||h || ff K (sl) , |||h||| T ) < CiV(||h || ff K (sl) , |||h||| T ). 



Now suppose that (0 is satisfied. If p" < 0, then 



1 



\\Hmu { si) + 2P(typ\mHt)\\H«- HS i)+ j o pWhww 

< C Sl N(\\h Q \\ H K {sl) , |||h||| T ) + C<$i|||h||lr- 
If logp' has small total variation, then for some 82 <C 1, 



liJK-^Bi) 



ds 



(78) 



\\0°gp'y\\ LH0iOO) 

Since py/p 7 \\h\\ H K-i(§i) < |||h||| T , 

p(s) 2 p'(s)||h(s)|| 2 ? K-i (s i ) ds 



poo 


P"(s) 


Jo 


p'(s) 



ds < 62- 



< 



°° \P"{*) 



p'(s) 



ds 



\W\\t < ^2|||h'" 2 



thus we obtain (|75|) again with 61 replaced by 61 +62- 

Remark 6.1. By the energy estimate (|78[) . we see that the quantity ||h||jyK-i(gi) 

decays at the rate — — . This decay rate is slower than the decay of ||h||#2.5(§i) if 
PVP' 

p grows algebraically. Thus, our approach of establishing decay in the lower-order 
norm via Duhamel's principle provides better decay than energy estimates alone. 

7. The stability of the Hele-Shaw flow with injection 

Combining (p?)) . (|rjg)l and ([75]), by choosing Si > small enough we conclude 
that 

\\\h\f T ^civaihoiiifK^jihiiy 



||h||| T < C 



2p 



hollffKCS 1 ) + ll|h||| T + |||h||| T 



(79) 



for some p £ N, provided that condition or ([7]), as well as the assumption 
(|29|) . are valid. The constant a is chosen sufficiently small so as to absorb certain 
error terms on the right-hand side of our energy estimates by the left-hand side 
energy terms. In other words, as long as all of the constants we use in our elliptic 
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estimates and the Sobolev embedding theorem are fixed, the maximum of a is a fixed 
computable number which can be chosen independent of ||ho || jye (S 1 ) - Inequality (|79|) 
then implies that there exists e > small enough such that if ||ho || jjKfs 1 ) < e > by 
the continuity (in time) of h, 

|||h||| T < a. (80) 

This suggests that as long as the solution exists, h has to satisfy the estimate (|5U|) . 
and this establishes Theorem 11.11 and 11.21 
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